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DESCRIPTION  OF  COMPUTER  PROGRAM 

Rotor  performance  and  loads  for  this  project  were  calculated 
by  Rotorcraft  Flight  Simulation  with  Aeroelastic  Rotor 
Representation,  Program  C81.  This  pregram  was  delivered  to 
the  U.  S.  i'rmy  under  Contract  DAAJ02-70-C-0063.  The  program 
is  completely  described  in  Reference  (1).  Following  is  a 
brief  summary  of  the  program. 

Program  C81  is  a general-purpose  rotorcraft  flight  simulation 
program  capable  of  analyzing  an  entire  rotorcraft  in  level 
or  maneuvering  flight.  However,  for  the  purposes  of  this 
study,  only  the  rotor  simulation  capability  was  used.  The 
rotor  simulation  is  divided  into  a trim  procedure  and  a loads 
calculation  portion.  The  trim  procedure  iterates  on  the 
independent  variables  (collective  pitch,  fore  and  aft  and 
lateral  cyclic  pitch,  and  fore  and  aft  and  lateral  lapping 
angles)  to  balance  the  forces  and  moments  acting  on  the  rotor. 
For  this  project,  the  program  trim  procedure  was  modified  to 
iterate  to  those  control  positions  necessary  to  give  zero  fore 
and  aft  and  lateral  flapping.  After  a trimmed  flight  condition 
is  obtained,  the  rotor  loads  calculations  are  performed  as 
described  below. 

The  rotor  loads  procedure  is  a fully  coupled  time-variant 
aeroelastic  analysis  based  on  the  medal  technique.  Modes 
and  natural  frequencies  must  be  calculated  in  a separate 
program  and  input  into  Program  C81.  The  modes  and  natural 
frequencies  used  in  this  study  were  generated  by  a program 
based  on  the  Myklestad  procedure. 

for  a specified  rotor  speed,  the  Ffyklestad  program  generates  a 
set  of  fully  coupled  modes  having  beam,  chord,  and  torsional 
components  as  a function  of  control  system  stiffness;  blade 
twist;  beam,  chord,  and  torsional  stiffness;  weight  distribu- 
tion; beam  and  chord  moments  of  inertia;  neutral  axis  and 
center-of-gravity  locations  at  twenty  radial  stations;  cen- 
trifugal stiffening;  centrifugal  restoring  torsional  moment; 
and  hub  impedance. 

Rotor  airloads  are  calculated  by  C81  at  20  radial  stations  for 
each  azimuth  location  as  a function  of  the  local  angle  of  at- 
tack and  local  Mach  number  to  include  the  effects  of  stall  and 
compressibility.  The  effects  of  both  rigid  body  and  elastic 
displacements  and  velocities  are  included  in  these  calculations. 
Also  included  at  the  user's  option  are  nonsteady  aerodynamic 
and  yawed  flow  effects.  Airfoil  lift,  drag,  and  pitching 
moment  coefficients  may  be  input  as  tables  of  discrete  values 
or  in  equation  form. 


InerL'.ia  forces  due  to  the  unsteady  motion  of  the  frame  of  ref- 
erence and  calculated  airloads  are  summed  to  form  the  forcing 
function  portion  of  the  modal  equations.  The  total  time- 
variant  response  of  the  blade  is  then  obtained  by  summing  the 
response  of  each  input  mode  shape.  Bending  moments  are  com- 
puted by  summing  the  products  at  each  radial  station  of  modal 
bending  moment  coefficients  and  modal  participation  factors. 

The  bending  moment  coefficients  for  each  mode  are  computed  as 
a function  of  blade  mass,  rotor  spaed,  mode  shape,  and  mode 
natural  frequency  by  considering  the  moments  generated  by  the 
inertia  forces  and  centrifugal  forces. 

Rotor  performance  and  blade  bending  moment  information  is  pre- 
sented at  the  end  of  each  Program  C81  simulation.  The  perform- 
ance information  is  given  in  both  dimensional  and  nondimension- 
al  form.  Bending  moments  are  presented  as  sine  and  cosine  har- 
monic analysis  components.  A complete  time  history  of  any 
rotor  performance  and  load  calculations  may  be  obtained  by 
entering  the  maneuver  portion  of  the  program.  This  section 
numerically  integrates  the  equations  of  motion  to  give  the  time 
history,  which  may  be  plotted  either  as  printer  plots  or  as 
CALCOMl?  plots. 


16 


H 


r»  • 


DESCRIPTION  OF  BLADE  PROPERTIES 

The  blade  dynamic  and  aerodynamic  data  were  obtained  by 
Sikorsky  Aircraft  Corporation  under  Contract  DAAJ02-72-C-0026. 
The  complete  description  of  the  experimental  techniques  is 
described  in  Reference  (2).  These  data  are  summarized  below. 

DYNAMIC  PROPERTIES 

Tlie  model  rotor  blades  examined  in  this  project  were  scaled 
from  the  Sikorsky  H-34  helicopter  rotor  system.  The  geometric 
scale  factor  applied  was  6.109.  The  four  blades  are  described 
as  follov7s: 


1. 

Fiberglass 

spar,  zero  cwist 

2. 

Fiberglass 

spar,  -8°  uniform  twist 

3. 

Fiberglass 
flap  (over 

spar,  zero  twist  with  5°  trailing-edge 
aft  207o  of  chord) 

4. 

Aluminum  spar,  zero  twist 

All  blade  sets  had  the  same  nominal  weight  and  inertia  distri** 
butions.  The  aluminum  blades,  however,  were  approximately 
three  times  stiffer  than  the  fiberglass  blades.  Blade  stiff- 
ness, inertia,  and  geometry  information  are  given  in  Tables  1 
through  4. 

The  tabulated  blade  stiffness  and  inertia  properties  were  ap- 
portioned in  twenty  equal  segments  to  match  the  input  format 
of  Program  C81.  These  equal  segment  length  data  are  shown  in 
Tables  5,  6,  and  7.  Calculated  frequencies  and  generalized 
inertias  are  given  in  Table  8.  Blade  inertia  and  center  of 
gravity  values  were  apportioned  on  a pro  rata  basis.  The 
unequal  segment  length  stiffness  values  were  combined  into 
equal  segment  length  values  by  using  the  following  relation: 


>■1*  ‘■2  + ■■■  * Vl 


El 


EQ 


(1) 


where 


-EQ 


an  equivalent  area  moment  of  inertia  for  a beam 
composed  of  n subsequents  having  different  area 


moments  of  inertia,  I. 
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I^,  and  lengths 
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This  procedure  yields  an  equivalent  beam  segment  which  rotates 
an  amount  equal  to  that  of  a beam  made  up  of  the  n subsegments 
when  both  beams  are  subjected  to  the  same  moment. 


AERODYNAMIC  PROPERTIES 


Two  sets  of  aerodynamic  data  are  presented  for  this  project. 
These  are  the  model  two-dimensional  airfoil  data  and  the  mod- 
ified airfoil  data  used  for  loads  predictions  (explained  and 
presented  in  the  Model  Rotor  Aerodynamic  Investigation  section 
of  this  report). 


The  model  lift,  drag,  and  pitching  moment  coefficients  are 
presented  in  Figures  1,  2,  and  3 for  the  symmetrical  (un- 
flapped) blade  section  and  in  Figures  4 5,  and  6 for  the 

trailing-edge  flap.  The  full-scale 
0012  litt  table  characteristics  used  in  the  modified  airfoil 
table  are  presented  in  the  Model  Rotor  Aerodynamic  Investiga- 
tion section  of  this  report. 


PRETEST  PREDICTIONS 
preliminary  CALCULATTONq 


M . -4U , M - U.438,  and  Ci /o  = 0.045  Theco 
points  \n: ire  selected  because  t?iey  appeared  to  be  reoresenta- 
tive  of  the  conditions  contained  in  Reference  (3). 


cuTacl™s'’’'°ThrrJl';ir^  "a  initial  performance  cal- 

cuiacions.  The  results  obtained  are  shown  in  Figure  7 whirh 

schedS^^d  C^/i°®liuefwLr«dSc:d 

line  conditions  were  then  Ci/o  = 0.049  and  cI/TL  0^1? 
scIib^OOO  calculations  made^uslng  the  full- 

?efiltf fL  “ =§1  nnd  meSsured 

resuitstor  the  fiberglass  blade  with  -8  twist.  These  twn 

Items  will  be  discussed  later  in  the  text. 


VARIATIONS  ABOUT  BASELINE  CONDITIONS 


One  hundred  sixty- two  C81  trim  cases  were  completed  to  match 

perfo"r^e1"ab1uf?h;T"1-^"\^  ^-e'clLl  werr'"" 

baseline  trim  points  for  the  fiberglass 
determine  the  sensitivity  of  pir? 
formance  and  loads  to  variations  (+2  ) in  mast  tilt  coller- 
^ve  pitch,  longitudinal  and  lateFal  cyclirpitch: 

^ese  cases  are  analyzed  in  the  Initial  Rotor^Performance 

Sberg?asr  baselines  f??  the 

ergiass,  -8  twist  blade  are  given  in  Figure  8. 


Additional  cases  were  calculated  for  each  baseline  condition 
in^he^o''®®  °ther  sets  of  blades  wi.rch  were  not  included 
in  Figures  aS'll. presented 
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CORRELATION  CRITERION 

A correlation  criterion  must  exist  to  judge  comparisons  of  ex- 
perimental data  and  calculated  predictions.  This  section  de- 
scribes the  criterion  used  to  quantify  the  results  of  this 
pro j ect . 

Performance  and  bending  moment  data  may  be  compared  point  to 
point  or  statistically,  using  curve  fits  to  relate  results. 
Each  method  has  its  benefits  and  drawbacks.  Point-to-point 
comparisons  yield  readily  understandable  results  for  each  set 
of  points  examined,  but  they  do  not  indicate  trends.  A sta- 
tistical representation  compacts  large  amounts  of  information 
and  describes  trends;  but,  when  taken  alone,  it  tends  to  hide 
individual  deviations.  Therefore,  use  of  both  procedures  is 
advisable. 


A quality  factor,  abbreviated  q,  is  used  with  the  point-to- 
point  and  statistical  comparisons.  For  point-to-point  com- 
parisons, a percentage  of  difference  is  used,  given  by 


q 


I measured  - computed 
measured 


X 100% 


(2) 


The  statistical  approach  involves  a different  definition, 
which  is  given  later  in  this  section.  Ranges  of  q are  given 
in  Table  9 • These  values  were  picked  from  experience  con- 
sidering the  ranges  of  accuracy  needed  by  the  designer. 

The  basis  of  the  statistical  method  is  to  plot  like  quantities 
(having  the  same  independent  parameters,  m-,  gQ,  ®3/a>  and 
a^i)  of  measured  data  versus  computed  predictions  and'^ro  fit 
a first-order  least- squared  error  criterion  curve  fit  through 
the  points.  The  correlation  criterion  may  then  be  applied  to 
the  parameters  describing  the  fitted  line.  The  line  has  three 
properties  which  describe  its  relation  to  the  data.  These  are 
its  slope,  intercept,  and  the  scatter  of  data  about  the  line. 
If  measured  data  (y)  is  plotted  against  computed  predictions 
(x),  a line  fitted  through  the  points  using  the  least- squared 
error  cri.terion  will  have  the  following  form: 

y'=ax+b  (3) 


where 


y'  = the  predicted  value  of  y 


a = 


2..XV  - nx 


Fx  - nx 


and 


b = l_  y^y 

Ix'^  - nx"^ 


n = the  number  of  samples 


A correlation  coefficient  may  be  defined  as 
S 2 

2 1 
r = 1 


V 


where 


r = the  linear  correlation  coefficient 

2 

= I ^ , the  variance  of  the  estimate 

e n-  i 


V = 


n 


= I 


the  mean  of  the  data 


, the  variance  of  the  data  about 
its  mean. 


(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 


Toe  quantity  may  be  interpreted  as  the  percentage  of  the 
variance  (S  2)  which  is  contributed  by  the  fitted  linear 
relation  fo¥  y' . It  is,  therefore,  a statistic  which  de- 
scribes the  usefulness  of  the  fitted  line  for  estimating  pur- 
poses. 

Perfect  correlation,  according  to  the  preceding  definitions, 
is  represented  by  a = 1,  b = 0,  and  r = 1.  The  terms  a and 
r are  independent  of  the  magnitude  of  the  variables;  however, 
b is  magnitude  dependent.  Therefore,  b is  norroalized  by  the 
average  value  of  the  independent  variable  and  subtracted  from 
one  so  that  unity  will  continue  to  represent  perfect  correla- 
tion. The  term  b is  replaced  by 


b'  = 1 - b/x 


(12) 


The  quality  factor  may  now  be  defined  (in  a percentage  of  error 
sense)  to  be 
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(l“a)  X 100%  for  a <•  1 


or 


% = (a-1)  X 100%  for  a > 1 


q^  - d-b')  X 100%  = b/x  X 100% 
qj-  = d - r)  X 100% 


S?t^r?“g%vfr[r?ab?;  r 


(13) 

(14) 

(15) 


"Introduction  to 
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INITIAL  ROTOR  PERFORMANCE  CORRELATION 


Because  of  the  large  quantity  of  data,  and  the  possibility 
that  the  actual  wind  tunnel  test  points  would  not  match  those 
used  in  the  previous  calculations,  statistical  techniques 
were  used  to  analyze  the  data.  The  following  data  were 
subjected  to  a multiple  correlation  analysis: 


Independent  Variables 
Vst’  ®3/4’ 


Dependent  Variables 


An  , Bi  , O-r 


C^/o,  Cy/a 
C;l/°,  C^O.  Cy= 


For  a specified  dependent  variable,  the  multiple  correlation 
analysis  picks,  in  order  of  descending  importance,  the 
independent  variables  that  may  be  used  in  a linear  combina- 
tion to  express  the  dependent  variable.  Within  the  162 
computer  cases  run,  the  rotor  lift  can  be  statistically 
represented  by 


Cf/o  = .00824  + .004  + .00693  6^/^ 


(16) 


The  multiple  correlation  coefficient  for  the  above  equation 
is  r = 0.951,  which  means  that  (.951)“^  x 100  = 90.57o  of  all 


variations  of  Cy /o  can  be  represented  by  variation  in 


and  0 


3/4* 


The  standard  error  is  S = 0.009,  which  can  be 

e 


interpreted  as  a confidence  limit  that  the  computed  Cj/o 


will  be  within  +0.009  of  the  value  predicted  by  the  above 
equation  in  65  percent  of  the  cases  and  within  +2x(0.009)  = 
+0.018  in  95  percent  of  the  cases. 


All  measured  data  for  the  fiberglass  -8  twist  blade 
were  also  subjected  to  the  multiple  correlation  analysis. 
Ttie  results  were 


n /o  = .01854  + .00456  <^rT,p-n  + .0066  6-3// 
L meas  -i-Fr  j/4 


(17) 


with  a correlation  coefficient  of  r = 0.911 


The  calculated  and  measured  lift  coefficients  may  be  compared 
by  examination  of  the  above  relations.  The  partial  deriva- 
tives giving  the  rate  of  change  of  lift  with  respect  to  tip 
path  plane  attack  angle  and  collective  setting  are  within 
12  percent.  The  difference  between  the  constant  terms  is 


in  the  2-D  wind  tunnel  data. 


due  to  the  low  values  of  C, 


max 


The  analysis  was  performed  on  the  calculated  and  measured 
power.  The  results  were  as  follows: 


Gr,/o  = 0.00006  + .00069  8./,  + .00010 

P'  meas  IPP 


S/"calc 


r = 0.928 
= 0.00020  + .00061  0 
r = 0.960 


3/4 


+ .00015  a 


TPP 


(18) 

(19) 


It  should  be  recalled  that  the  correlation  analysis  is  based 
on  linear  theory  and,  therefore,  any  nonlinear  relationship 
among  the  variables  is  deleted.  Since  the  multiple  correla- 
tion analysis  is  free  to  select  any  of  nine  independent 

is  significant  that  collective  pitch  and  tip 
attack  were  always  selected  as  the  best 
rotor  thrust  and  power.  The  fact  that  these 
made  for  both  the  measured  and  calculated 
validate  the  mathematical  model  of  the  rotor. 


variables,  it 
path  angle  of 
predictors  of 
selections  were 
values  tends  to 


The  statistical  representation  of  the  measured  and  calculated 
gross  rotor  performance  is  shown  in  Figures  11  and  12.  The 
computed  results  used  in  this  analysis  are  tabulated  in  Table 
9.  The  computer  subroutines  used  in  the  multiple  regression 
analysis  are  described  in  Reference  (5). 
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SENSITIVITY  ANALYSIS 


An.  analysis  was  developed  to  relate  the  sensitivity  of  pre- 
dicted results  to  input  parameter  variations.  Twenty- six 
individual  variations  of  significant  input  parameters  were 
made  for  the  two  baseline  correlations  (Cl/o  = 0.049  and 
Cl/o  = 0.027)  using  the  model  aerodynamic  data.  However, 
as  explained  in  the  Performance  Comparisons  Section  of  this 
report,  the  program  could  not  accurately  represent  the  model 
rotor  tunnel  behavior  using  the  model  steady  aerodynamic  data. 
This  fact,  coupled  with  apparent  discrepancies  in  the  model 
two-dimensional  data  apart  from  any  program  considerations, 
made  the  sensitivity  analysis  calculations  generally  suspect. 
Consequently,  these  calculations  are  not  presented  in  this 
report.  However,  the  mathematical  procedure  for  the  sensi- 
tivity study  is  considered  valid  and  very  useful  for  future 
studies.  Iherefore,  the  analysis  is  presented  here  for  in- 
formation only. 


If  the  standard  deviation  for  an  input,  x-[,  ip  o^,  then  the 
standard  deviation  of  an  output  quantity,  y,  may  be  written 


^ i=l  V 


(20) 


where 


n = number  of  input  variables,  x 


(21) 


^ = rate  of  change  output,  y,  with  respect 

^i  to  input,  Xj^. 


(22) 


To  employ  this  method,  a value  must  be  assigned  to  each  input 
parameter  standard  deviation.  These  values  should  be  assigned 
from  the  results  of  a large  sample  test  of  the  values  of  each 
quantity.  Lacking  this  information,  a judgment  must  be  made 
to  relate  the  standard  deviation  to  confidence  levels  (given 
in  plus  and  minus  percent  error  about  a nominal  value). 


If,  after  the  determination  of  oy,  the  computed  value  and 
measured  value  of  y are  compared  and  are  such  that 


y - y, 


measured 


(23) 


then  there  is  a 95  percent  probability  that  the  variation 
between  the  measured  and  computed  values  is  due  to  input 
parameter  variation. 


The  method  given  here  follows  Bendat  and  Piersol  (Reference  6). 
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POSTTEST  PERFORMANCE  CORRELATION 


Perfordiance  comparisons  of  calculated  and  measured  results 
were  made  for  the  following  cases: 


Fiberglass 

blade. 

-8" 

twist,  6.^  = 0^^ 
r 

i-i 

0.299 

^1,90 

0.408 

0.400 

0.435 

0.502 

0.467 

Fiberglass 

blade, 

0^' 

twist,  t-p  = 0° 

0.299 

0.408 

Fiberglass 

blade , 

0‘^ 

twist,  = 5^ 

0.299 

0.40S 

Comparisons  of  measured  and  calculated  results  are  shown  in 
Figures  14  through  33.  All  calculated  results  sho\^7n  in  these 
figures  were  made  using  the  model  airfoil  data.  Performance 
comparisons  for  the  aluminum,  0°  twist  blade  are  not  shown 
since  the  analysis  used  in  this  portion  of  the  study  is  a 
quasi-static  analysis  (using  steady  and  1/rev  blade  response) 
which  does  not  yield  results  dif f arent  from  those  for  the 
fiberglass,  0 twist  blade.  Basic  comparisons  are  presented 
as  plots  of  lift  versus  control  plane  angle  of  attack,  power 
versus  lift,  measured  lift  versus  calculated  lift,  and  measur- 
ed power  versus  calculated  power.  Measured  lift  versus  cal- 
culated lift  and  measured  power  versus  calculated  power  are 
not  shown  for  the  fiberglass,  -8^  tvjist  blade  at  the  k = 

0.502,  Mj^  gQ  = 0.467  condition  because  insufficient  measured 
data  were ’ available  for  a statistical  analysis. 

Examination  of  the  plots  shows  that  measured  and  calculated 
lift  values  are  in  close  agreement  at  low  collective  values. 
Calculated  power  is  generally  lower  than  measured.  As  collec- 
tive pitch  increases,  the  calculated  lift  shows  stall  char- 
acteristics at  a much  lower  pitch  angle  than  the  measured 
lift.  These  stall  characteristics  are  best  compared  in  the 
power  versus  lift  plots.  These  figures  show  calculated  stall 
to  be  extremely  abrupt.  For  instance,  in  Figure  15  a maxi- 
mum calculated  lift  coefficient  of  Cj^/o  = 0.07  has  been  reach- 
ed while  calculated  power  increases  asymptotically.  Calcula- 
ted stall  is  most  severe  for  the  blade  set  with  the  5 trail- 
ing-edge  flap  (Figure  31),  as  expected.  The  measured  stall 
characteristic  for  this  case  is  hardly  in  evidence,  no  more  so 
than  for  the  unflapped  blade  at  the  same  condition  (Figure  27). 
Figures  18  and  19  illustrate  the  effect  of  the  Program  C81  un- 
steady aerodynamic  and  yawed  flow  representations  on  the  cal- 
culated lift  and  power  for  the  fiberglass  blade  with  -8° 
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pimpp 


Figures  14  and  15).  Though  t.hese  tlftct** 
incrt.a^e  the  calculated  lilt  at  high  collective  settings 

Tabli  abruptly  than  in  the  measured  rfsilL 

SuuliLS  fflf  of  measured  Ind 

ca-icuiated  iitt  and  power  and  also  grades  the  correlation 

according  to  the  correlation  criteria  given  in  Table  9. 

rnvLS^  ^I'^-lated  reLlts  datar^^urthT^'' 

fo^^  indicated.  This  investigation  took  the 

form  of  a brief  correlation  study  for  a full-scale  rotor 

rotors  ^^The  the  aerodynamics  of  model 

two  sections  studies  are  given  in  the  next 


i 


27 


FULL-SCALE  CO^^PAJ  .SONS 

Program  C81  was  used  to  predict  performance  and  bending  moments 
for  a full-scale  H-34  rotor  system.  The  experimental  data  for 
this  comparison  are  published  in  USAAVLABS  Technical  Report 
68-3  (Reference  3)  and  NASA  TND-4632  (Reference  7).  All  blade 
dynamic  properties  for  these  comparisons  were  estimated  from 
the  figures  and  plots  given  in  TR68-3.  The  aerodynamic  prop- 
erties were  those  of  standard  0012  tables  used  at  Bell  Heli- 
copter. 

Performance  comparisons  for  various  mast  tilt  angles  are  given 
in  Figures  34  through  39.  Figures  34  through  36  compare  meas- 
ured and  calculated  lift,  drag,  and  horsepov;er  with  collective 
pitch  as  the  independent  parameter,  while  Figure  37  shows  meas- 
ured and  calculated  power  versus  lift.  As  may  be  seen,  the 
calculated  performance  parr  -iters  show  very  good  qualitative 
agreement  with  the  measured  data.  Lift  error  is  a maximum  of 
900  lb  for  the  0°  mast  tilt  case  at  ^3/4  = 2 , amounting  to  18 
percent  of  the  measured,  and  decreases  with  increasing  collec- 
tive pitch.  Drag  error  is  about  500  lb  or  approximately  30 
percent  of  measured  for  the  10°  mast  tilt  case.  Maximum  power 
error  is  approximately  200  hp  or  about  16  percent  of  measured. 
Figure  37  illustrates  that  the  full-scale  rotor  stall  char- 
acteristics are  predicted  very  accurately. 

Figures  38  and  39  compare  experimental  lift  to  calculated  lift 
and  experimental  power  to  calculated  power.  This  method  of 
presentation  allows  application  of  the  correlation  method  de- 
scribed previously.  First-order  least-squared  error  criterion 
prediction  lines  for  each  of  these  plots  are  described  by 


= 0-0118  + 0.923  C,/o  . 
L'  meas  L'  calc 


(24) 


with  r = 0.989  and  = 0.061 


o = 0.0000583  + 1.014  Qjo  , 

P meas  P'  calc 


(25) 


with  r = 0.957  and  Cp/o  = 0.00239. 

The  quantities  a(a'),  b(b’),  and  r (slope,  normalized  inter- 
cept, and  cor'reiation  coefficient  as  defined  previously)  for 
lift  and  power  may  be  expressed  as 


Lif  t : 


a = 0.923  (q  = 7.8%) 


(26) 


Power: 


' ‘(v“) 


Q.0118\ 

J7a6i“j 


= 0.807 


(27) 


(q  = 19.37c) 
r = 0.989  (q  = 1.17,) 
iv: 

a = 1.014  (q  = 1.47o) 

(q  = 2.447c) 
r = 0.957  (q  = 4.3‘/c) 


0.0000583 


■)  = 0.97t 


(28) 


(29) 

(30) 


(31) 


Application  of  the  correlation  criterion  stated  in  Table  9 
grades  the  lift  slope  correlation  as  good,  lift  intercept  as 
poor,  and  lift  correlation  coefficient  as  excellent.  Appli- 
cation of  the  criterion  to  power  grades  the  power  slope  as 
excellent,  power  intercept  as  excellent,  and  power  correlation 
as  excellent. 

Figures  40  and  41  compare  measured  and  computed  bending  stresses. 
The  Program  C81  beamwise  stress  predictions  are  higher  in  the 
midspan  region  than  the  measured,  but  are  in  good  agreement  in- 
board of  r/R  = 0.35  and  outboard  of  r/R  = 0.70.  However,  the 
predicted  chordwise  stresses  are  low  in  the  inboard  regions  and 
quite  good  outboard  of  approximately  r/R  = 0.40. 


MODEL  ROTOR  AERODYNAMIC  INVESTIGATION 


This  section  presents  a brief  examination  of  the  performance 
correlation  problems  described  in  the  Posttest  Performance 
Correlation  section.  The  model  rotor  aerodynamics  were  exam- 
ined and  a modified  airfoil  developed  which  produced  reason- 
able performance  correlation. 


Rotor  stall  character.i.stics  are  compared  in  Figures  42  and  43. 
Figure  42  presents  measured  power,  computed  power  (using  model 
data  tables),  and  measured  oscillatory  moment  (at  r/R  = 0.20) 
plotted  against  lift.  The  case  presented  is  at  a = 2°  to 
represent  a forward  flight  condition.  The  plotted  power  fig- 
ures are  modified  by  removal  of  induced  power  and  parasite 
power  to  isolate  aerodynamic  effects.  These  terms  are  approx- 
imated by  simple  momentum  theory  and  are  not  exact.  They  are 


Induced  power  coefiicient  = Cpjj^/ci  = (Cj^/o)2  o/2n 


Parasite  power  coefficient  = Cp  /°  = m(C  /o) 

P ^ 


(where  C /a  is  the  rotor  x-force  coefficient  in  the  wind  axes). 

X 


Figure  42  clearly  illustrates  that  measured  power  and  measured 
torsional  moment  "break"  together,  further  evidence  that  the 
phenomenon  observed  is  stall.  The  computed  power  breaks  well 
before  the  measured  power.  This  observation  is  considered 
further  evidence  that  computations  made  with  the  model  airfoil 
tables  do  not  represent  the  rotor  as  it  behaves  in  the  tunnel. 
Figure  43  is  another  power  versus  lift  plot  which  presents 
comparisons  of  measured  data  and  predictions  computed  using 
Program  C81  and  Program  F-35  (Ref.  8),  a Bell  Helicopter  rotor 
performance  program.  Both  computations  were  made  using  the 
model  airfoil  data  tables  to  check  the  programs  against  each 
other.  As  may  be  seen,  both  programs  compute  essentially  the 
same  stall  characteristics. 


The  results  described  above,  the  stall  characteristics  observed 
in  the  Posttest  Performance  Correlation,  and  the  full-scale 
comparisons  all  suggest  that  the  performance  correlation  prob- 
lems were  caused  by  the  stall  characteristics  of  the  model  air- 
foil data.  It  is  well  known  that  the  Reynolds  number  and  the 
Mach  number  can  have  significant  effects  on  blade  aerodynamic 
coefficients.  The  2-D  airfoil  tests  performed  by  Sikorsky  were 
designed  specifically  to  obtain  the  necessary  steady-state  data 
to  be  used  in  the  flight  simulation  program.  These  data  are 
presented  in  Figures  1,  2,  and  3.  Of  particular  interest  is 
the  maximum  lift  coefficient  as  a function  of  Mach  number 

shown  in  Figure  44.  Also  shown  ror  comparison  are  values  of 


Ct^  from  synthesized  airfoil  tables  developed  from  a Bell 
in^ouse  hovering  rotor  test  and  a Sikorsky  hovering  rotor 
test  (Reference  9).  The  Bell  test  was  conducted  with  a two- 
bladed  2.25  in.  chord,  14  in.  radius  blade.  Three  twist  com- 
binations and  three  tip  Mach  numbers  were  examined  and  two- 
dimensional  data  were  developed  that  correlated  the  measured 
performance  for  all  three  rotors  at  all  three  Mach  numbers. 

This  approach  has  been  used  by  NACA  (Reference  10)  to  synthe- 
size airfoil  data  from  test  results.  Both  sets  of  synthesized 
data  contain  significantly  higher  maximum  lift  coefficient 
values  when  compared  to  the  model  two-dimensional  data.  The 
conclusion  must  be  that  the  model  data  is  unreasonably  low  or 
that  unusual  effects  caused  the  synthesized  data  to  be  unex- 
pectedly high. 

Ttie  model  two-dimensional  airfoil  data  was  used  to  predict 
performance  at  p = 0.299,  Mi  90  ~ 0.408  and  = 0°  as  shown 
in  Figure  7.  When  the  model ’data  are  used,  stall  occurs  at 
approximately  Cl/o  = 0.075.  These  values  are  considerably  be- 
low the  measured  lift  values  shown  for  this  condition.  If 
full-scale  0012  airfoil  data  are  used,  the  stall  point  is 
increased  to  0.106.  Comparison  between  the  measured  and  full- 
scale  lift  values  indicates  that  the  full-scale  data  more 
nearly  approximates  the  behavior  of  the  rotor  in  the  tunnel. 

The  exact  reason  for  the  discrepancy  between  measured  and  cal- 
culated results  is  not  well  understood;  however,  it  may  be 
related  to  "centrifugal  pumping".  Centrifugal  pumping  refers 
to  the  tendency  for  the  centrifugal  force  to  push  the  turbulent 
boundary  layer,  which  occurs  with  blade  stall,  out  of  the 
blade.  As  the  turbulent  volume  is  removed,  the  flow  reattaches 
to  the  airfoil  and  thus  the  stall  conditions  are  removed.  The 
effect  may  be  more  significant  for  model  rotors  than  for  full- 
scale  rotors  because  of  the  increased  centrifugal  force  field. 

Since  the  model  behavior  in  the  tunnel  seemed  to  be  best  de- 
scribed by  full-scale  rotor  aerodynamic  characteristics,  a 
modified  airfoil  table  was  assembled  to  permit  bending  moment 
calculations  to  be  made.  The  table  used  for  the  three  blades 
without  the  trailing  edge  flap  contained  a full-scale 
(Rn  = 12  X 10^)  lift  table,  the  model  drag  table  with  an  incre- 
ment added  uniformly  to  all  values  (4Cp  = 0.003),  and  a null 
pitching  moment  table.  Figure  45  presents  the  lift  character- 
istics of  the  full-scale  0012  tables  used  in  the  modified 
tables.  Drag  characteristics  may  be  obtained  from  the  model 
drag  properties  shown  in  Figure  2 by  adding  the  increment. 
Figures  46  through  50  present  representative  effects  of  the 
modified  airfoil  on  performance  calculations.  Different 
increments  were  used  for  the  flapped  airfoil.  These  were 
4Cl  = 0.15  added  to  the  full-scale  lift  table  and  4Cp  = 0.005 
added  to  the  model  drag  table. 
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Applying  the  standards  presented  in  Table  9,  the  lift  corre- 
lation i?  graded  as  good,  excellent,  and  excellent  for  its 
slope,  intercept,  and  correlation  coefficient.  Sirr.ilarly, 
power  correlation  is  graded  as  fair,  good,  and  excellent  on 
its  slope,  intercept,  and  correlation  coefficient. 


The  modified  airfoil  was  used  because  it  allowed  the  calcula- 
tion of  reasonable  performance  parameters.  A general  technique 
for  adjusting  two-dimensional  model  airfoil  data  to  represent 
actual  operating  airfoil  characteristics  has  yet  to  be  developed. 
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POSTTEST  ROTOR  BENDING  MOMENT  CORRELATION 


Load  comparisons  were  made  by  computing  blade  response  at  cer- 
tain specified  conditions  for  all  four  blade  sets.  The  condi- 
tions were  selected  to  present  a range  of  lift  and  p,  qq 

combinations.  The  cases  examined  were  designated  in  the  Sikor- 
sky wind-tunnel  test  as  Conditions  25,  44,  and  68.  Because 
each  condition  number  represents  slightly  different  operating 
conditions  for  each  blade,  the  values  of  lift  and  power  coeffi- 
cients corresponding  to  each  blade  set  are  listed  in  Table  11. 
This  table  also  shows  the  computed  lift  and  power  for  each  con- 
dition and  blade. 

Load  computations  were  made  using  the  modified  airfoil  table 
discussed  in  the  Model  Rotor  Aerodynamic  Investigation  section 
of  this  report.  For  each  blade  set  and  condition,  the  computed 
lift  and  first  harmonic  flapping  components  were  brought  as 
close  as  practicable  to  the  measured  lift  and  flapping.  A col- 
lective sweep  using  a quasi-static  analysis  (steady  and  first 
harmonic  blade  response  only)  was  made  to  obtain  the  value  of 
lift  required  by  the  target  condition.  The  Program  C81  quasi- 
static procedure  will  automatically  iterate  to  zero  first  har- 
monic flapping  component  values.  However,  when  a fully  elastic 
computation  is  made,  the  values  of  lift  and  flapping  angles  com- 
puted in  the  quasi-static  analysis  are  changed  slightly 'due  to 
aeroelastic  effects  (principally  changes  in  blade  angle  of 
attack  due  to  local  blade  velocity).  These  changes  required  a 
succession  of  simulations  to  obtain  the  desired  lift  and  flap- 
ping angles.  Consequently,  the  computed  values  do  not  exactly 
match  the  measured.  Also,  because  the  modified  airfoil  char- 
acteristics do  not  represent  the  model  behavior  adequately  in 
all  regimes,  there  is  frequently  a mismatch  between  computed 
and  measured  power. 

A choice  of  rotor  blade  elastic  modes  was  necessitated  by  a 
Program  C81  input  limit  of  six  modes.  Properties  of  the  modes 
used  are  shown  in  Table  8.  Computations  for  the  fiberglass 
blades  were  made  using  one  rigid-body  beamwise  mode,  one  rigid- 
body  chordwise  mode,  two  elastic  beamwise  modes,  and  two  elas- 
tic chordwise  modes.  For  these  blades,  the  selection  of  the 
second  elastic  chordwise  mode  at  approximately  7.5/rev  pre- 
cluded the  use  of  an  elastic  beamwise  mode  existing  at  approx- 
imately 6.7/rev.  This  choice  was  made  to  provide  a possibility 
of  higher  frequency  chordwise  response.  Computations  for  the 
aluminum  blades  were  made  using  the  first  existing  six  modes. 
Torsion  modes  were  omitted  for  ail  calculations  because  the 
modified  airfoil  had  no  pitching  moment  table.  Modes  were 
calculated  for  a rotor  speed  of  730  rpm  and  a root  collective 
pitch  of  12°.  The  boundary  conditions  applied  for  mode  shape 
calculations  were  those  of  a pinned  out-of-plane  hinge  and  a 
pinned  in-plane  hinge  located  at  r * 3 inches.  The  effect  of 
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the  inplane  lead-lag  damper  v;as  represented  in  Program  C81  by 
a damping  force,  proportional  to  inplane  velocity,  applied  at 
r/R  = 0.10.  Based  on  the  damper  rate  of  17  in. -Ib-sec/rad , 
the  Program  C81  damper  input  was  37.5  Ib-sec/ft. 


Measured  and  calculated  beam  and  chord  bending  moments  are  com- 
pared by  time  histories  and  by  harmonics  at  five  radial  sta- 
tions (r/R  = 0.20,  r/R  = 0.35,  r/R  = 0.45,  r/R  = 0.65,  and 
r/R  = 0.80).  Figures  51  through  74  present  the  time  histories 
for  the  twelve  cases  examined.  These  are  followed  by  Figures 
75  through  98  which  present  the  harmonic  analyses  of  the  same 
cases.  The  harmonics  are  given  in  magnitude  and  phase  form  as 
explained  below.  Tables  12  through  35  present  the  measured 
and  calculated  bending  moments  as  sine  and  cosine  harmonic 
components. 


Measured  and  calculated  load  harmonics  are  presented  in  Figures 
75  through  98  as  magnitudes  and  phases  plotted  against  radial 
station.  A load,  M,  may  be  represented  as  Fourier  series  sine 
and  cosine  components  in  the  form 


M(nC^t,r)  sin  ni<t  + B cos  nQt) 

n=0 


(32) 


where 


M = load  magnitude 
n = harmonic  number 
0 = rotor  speed 
t = time 

r = blade  radial  station 

A^  = sine  component  of  the  load  at 

radial  station  r 

B_  = B (r),  the  cosine  component  of  the  load 
at  radial  station  r. 


n 


The  above  relation  may  also  be  written  as  a magnitude  and 
phase  in  the  form 


M(nfit,r)  ='^  cos  (nfit  - 0^) 
n=0 


(33) 


where 


= A^^  + B.  ^ 


n 


n 
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Of  the  twelve  cases  examined,  calculations  using  the  -8°  twist 
blade  produced  the  best  results.  Examination  of  Figures  57 
and  58  shows  oscillatory  magnitudes  to  be  quite  good  for 
Condition  25.  The  beam  traces  for  this  condition  show  good 
magnitude  and  phase  relationships  for  the  entire  span.  Chord- 
wise  loads  are  good  from  r/R  = 0.35  to  r/R  = 0.80.  However, 
there  is  a decided  difference  in  magnitude  at  the  20-percent 
radius.  This  effect  is  very  evident  in  the  plots  of  the  chord 
load  harmonics  (Figure  82),  where  good  correlation  is  seen 
for  1/rev  loads  from  r/R  = 0.35  to  r/R  = 0.80,  and  for  the  2/rev 
and  3/rev  chord  load  components  for  the  entire  span.  The  most 
noticeable  difference  is  seen  in  the  1/rev  chord  loads  at 
r/R  = C.20,  where  the  calculated  load  is  approximately  50 
percent  of  the  measured.  Calculated  beam  and  chord  load  magni- 
tudes for  Condition  44  are  low.  Pnase  relationships  are  good 
for  the  1/rev  beam  bending  moments  over  the  entire  span  and 
improve  for  the  2/rev  loads  toward  the  tip.  Once  again,  a 
large  jump  in  chord  1/rev  load  is  observed  at  the  20-percent 
radius.  Condition  44  measured  chordwise  moments  show  a very 
high  4/rev  harmonic  content  which  is  not  calculated  and  also 
not  seen  in  either  Condition  25  or  Condition  68  for  this 
blade.  Oscillatory  magnitudes  for  Condition  68  are  low  over 
the  entire  span  for  both  beam  and  chord  loads.  Calculated 
beam  harmonics  are  low  at  all  frequencies.  Calculated  chord 
1/rev  magnitude  again  shows  a very  large  discrepancy  at  the 
20  percent  radius  but  shows  continuing  improvement  in  magni- 
tude and  phase  with  blade  radial  station.  Calculated  chord 
2/rev  load  is  low  at  the  20%  radius,  but  shows  almost  perfect 
m'lonitude  correlation  from  r/R  = 0.45  to  the  blade  tip. 
Calculated  chord  3/rev  and  4/rev  magnitudes  are  low  over  the 
entire  span. 

Calculated  loaus  for  the  fiberglass  blade  with  0°  twist  are 
lower  than  measured  for  all  three  conditions.  The  measured 
loads  also  exhibit  high  harmonic  content  not  seen  in  the 
calculated  loads.  The  measured  chordwise  loads  for  Conditions 
25  and  68  show  large  excursions  in  the  fourth  quadrant  of  the 
rotor  azimuth  which  are  indicative  of  stall.  These  are  not 
accompanied  by  a calculated  stall  prediction. 

Calculated  loads  for  the  remaining  two-blade  sets,  the  fiber- 
glass blade  with  the  trailing-edge  flap  and  the  untwisted 
aluminum  blade,  are  considerably  lower  than  the  measured  loads. 
Generally  calculated  loads  do  not  show  the  higher  harmonic 
content  of  the  measured  loads,  and  calculated  stall  indica- 
tions are  absent. 


Calculated  loads  are  characterized  by  the  absence  of  higher  har- 
monic content  and,  excepting  the  fiberglass,  -8  twist  blade  at 
Condition  25,  generally  low  magnitudes.  Measured  1/rev  chord- 
wise  loads  seem  to  be  always  increasing  in  the  inboard  direc- 
tion, toward  the  hinge.  This  tendency  would  seem  to  indicate 
some  moment  carryover  at  the  lead-lag  hinge.  The  high  harmonic 
content  of  the  measured  loads  may  indicate  the  presence  of 
blade  interaction  effects  which  are  not  presently  modeled  in 
Program  C81. 

EFFECT  OF  CYCLIC  RESOLUTION 

Another  factor  which  may  influence  the  bending  moment  correla- 
tion is  the  effects  of  the  cyclic  pitch  in  the  resolution  angle 
between  the  shaft  reference  system  and  the  blade  reference  sys- 
tem. Referring  to  the  description  of  the  mathematical  model, 
it  can  be  seen  that  the  blade  bending  moments  are  formed  by  sum- 
ming the  bending  moments  produced  by  each  mode.  At  the  begin- 
ning of  each  computer  run,  the  input  mode  shapes  are  used  to 
calculate  the  in-plane  and  out-of-plane  bending  moment  coeffi- 
cients for  each  mode.  During  the  time  variant  aeroelastic 
rotor  analysis,  the  response  of  each  mode  or  the  participation 
factor  is  obtained  from  the  numerical  integration  of  the  equa- 
tions of  motion. 

These  participation  factors  are  multiplied  by  the  previously 
mentioned  in-plane  and  out-of-plane  bending  moments.  The  bend- 
ing moments  are  then  resolved  through  the  steady  geometric 
pitch  angle  (collective  pitch  plus  structural  twist)  to  obtain 
the  beam  and  chord  bending  moments  which  have  been  presented  in 
this  report.  Thus  the  cyclic  pitch  is  not  included  in  the  res- 
olution angle.  Later  versions  of  C81  have  been  modified  to  in- 
clude the  cycl.ic  pitch  angle  in  the  resolution  angle. 

Figures  99  and  100  illustrate  the  effect  of  the  cyclic  resolu- 
tion angle  on  loads.  The  calculations  shown  in  these  figures 
were  made  using  a version  of  Program  C81  which  incorporates  the 
cyclic  pitch  angle  in  the  angle  of  resolution.  The  results 
shown  are  typical  for  beam,  and  chord  bending  moments.  The  most 
noticeable  effect  is  a general  increase  of  beamwise  moments  in 
the  midspan  and  outer  regions  of  the  blade. 

EFFECT  OF  UNSTEADY  AERODYNAMICS  AND  YAWED  FLOW 

The  Program  C81  unsteady  aerodynamics  were  employed  to  compute 
the  bending  moments  presented  in  Figures  101  and  102  for  the 
fiberglass  blade  with  0°  twist  at  Condition  25.  These  figures 
should  be  compared  with  Figures  51  and  52,  which  present  com- 
puted results  without  the  unsteady  aerodynamics.  Load  magni- 
tudes are  not  affected  significantly  in  either  the  beam  ov 
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chord  moments.  The  effect  does  produce  a damping  effect  which 
is  seen  as  a lagging  phase  shift  (moving  the  calculated  loads 
toward  the  measured;. 


Unsteady  aerodynamics  and  the  yawed  flow  simulation  were  used 
to  produce  the  representative  time  history  trace  shown  in 
Figure  103.  The  effects  do  not  modify  the  load  amplitudes  ap- 
preciably. However,  the  introduction  of  the  yawed  flow  repre- 
sentation does  tend  to  shift  the  calculated  phase  in  the  lead- 
ing direction,  cancelling  the  phase  shift  of  the  unsteady 
aerodynamics  effect. 


EFFECT  OF  TORSIONAL  MOMENT  CALCULATIONS 


The  effect  of  incorporating  a torsional  mode  shape  and  a pitch- 
ing moment  table  is  presented  in  Figure  104  for  the  fiberglass 
blade  with  -8°  of  twist  at  Condition  44.  The  inclusion  of  the 
7.36/rev  torsional  mode  required  the  removal  of  a chordwise 
mode  at  7.56/rev  due  to  the  Prograia  C81  input  limit  of  six 
modes.  The  pitching  moment  table  used  v;as  the  model  0012  table. 
Figure  104  shows  that  tnis  effect  does  not  appreciably  affect 
results.  This  figure  is  representative  of  the  order-of- 
magnitude  difference  seen  at  all  radial  stations. 


EFFECT  OF  IN-PLANE  SPRING 


Since  the  calculated  1/rev  chordwise  moments  at  r/R  = 0,20  were 
consistently  lower  than  the  measured  values,  the  possibility  of 
excessive  moment  carryover  at  the  lead-lag  hinge  was  considered. 
The  Sikorsky  data  package  for  this  program  states  that  possibly 
15  percent  of  the  in-plane  lead-lag  damper  value  could  be  con- 
sidered as  a spring  component.  Assuming  1/rev  as  the  basis  for 
calculating  this  spring  component  gives  a value  of  approximately 
0.3  ft-lb/deg.  Modes  and  natural  frequencies  incorporating 
this  in-plane  spring  value  were  used  to  compute  the  representa- 
tive time  history  trace  shown  in  Figure  105.  The  effect  of  the 
spring  is  to  shift  the  mean  value  of  the  moments  to  a slightly 
higher  level.  The  higher  1/rev  chordwise  moments  may  also  be 
due  to  moment  carryover  by  the  lead- lag  hinge  bearings. 


AERODYNAMIC  REPRESENTATION 


Correlation  problems  encountered  are  believed  to  be  due  mainly 
to  an  inadequate  aerodynamic  representation.  Tne  modified  air- 
foil used  for  load  predictions  apparently  did  not  represent  the 
airfoil  well  in  lift  and  drag  and  contained  no  pitching  moment 
table  since  there  was  no  method  for  determining  moment  coef- 
ficient values. 


MEASURED  DATA 

Some  problems  were  observed  in  the  data  used  for  comparison  in 
this  project.  Figures  106  and  107  show  representative  time 
histories  of  measured  model  beamwise  and  chordwise  moments. 

The  chordwise  moments  contain  strong  noninteger  per/rev 
components  which  may  affect  comparisons.  The  beamwise  moment 
time  history  is  more  consistent  but  still  contains  components 
which  may  be  significant  at  high  frequencies.  Though  non- 
harmonic response  may  occur  naturally  due  to  various  causes, 
this  behavior  is  believed  to  be  primarily  related  to  inplane 
moments  generated  by  the  rotor  drive  system.  Though  the 
time  histories  presented  herein  are  averages  of  measurements 
made  for  40  revolutions,  the  illustrated  behavior  may  still 
significantly  affect  results.  The  nature  of  the  data  indicates 
that  measured  results  should  be  analyzed  to  detect  and 
identify  noninteger  per/rev  responses. 
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CONCLUSIONS 


The  Rotorcraft  Flight  Simulation  Program  (C81)  can 
predict  performance  and  bending  moments  for  full-scale 
rotors . 

Since  the  Rotorcraft  Flight  Simulation  has  been  shown  to 
correlate  well  for  a similar,  but  full-scale,  rotor 
system,  the  problems  experienced  are  likely  due  to 
unsteady  and/or  Reynolds  number  and/or  other  scale 
effects  which  are  not  accounted  for  in  the  program. 

The  measured  model  rotor  performance  is  not  compatible 
with  the  model  two-dimensional  aerodynamic  coefficients. 

Mathematical  modeling  techniques,  derived  for  full-scale 
rotor  aerodynamics,  may  not  be  valid  for  model  rotors. 

The  sensitivity  analysis  which  was  performed  is  not 
valid  due  to  limitations  in  the  input  data  and/or  the 
computer  program. 

Simple  harmonic  analysis  of  measured  data  may  lead  to 
erroneous  conclusions  because  of  its  inability  to 
detect  the  noninteger,  or  transient,  components. 

The  correlation  criterion  developed  for  rotor  perform- 
ance behavior  is  not  applicable  to  model  rotor  bending 
moment  correlation. 

The  difficulties  associated  with  model  rotors  should  not 
preclude  the  careful  use  and  orderly  development  of  the 
analysis . 


RECOMMENDATION.q 
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Model  Aerodynamic  Lift  Coefficient  vs, 
Angle  of  Attack. 
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Pitching  Moment  Coefficient 


Figure  3.  Mode^AirfoU^atcMng  Moment  Coeffld, 


Angle  of  Attack,  or,  deg 


Figure  A.  Model  Aerodynamic  Lift  Coefficient  vs 
Angle  of  Attack  for  Flapped  Airfoil. 
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Model  Aerodynamic  Pitching  Moment  Coefficient 
vs.  Angle  of  Attack  for  Flapped  Airfoil. 
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Figure  8.  Baseline  Case  Fiberglass  Blade  with  -8  Twist 

M = 0.398,  gQ  = 0.440,  Z^/o  = 0.027. 
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Figure  8.  Concluded:  |a  = 0.298,  M-,  = 0.408,  C, /o  = 0.049 
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Figure  9.  Concluded:  |a  = 0.298,  M-,  qp.  = 0.408,  C,/a  = 0.049 
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Figure  11.  Baseline  Case  Aluminum  Blade  with  0°  Twist 
M = 0.398,  M,  QQ  = 0.440,  C, /o  = 0.027. 
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Figure  11.  Concluded;  m = 0.298,  M-,  qp>  = 0.408,  C-, /o  = 0.049 
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Figure  13.  Comparison  of  Regression  Analysis  Results 
for  Measured  and  Calculated  Power. 
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Figure  17.  Experimental  Power  Coefficient  vs.  Calculated  Power 
Coefficient,  Fiberglass  Blade,  -8°  Twist,  ^l=0.299, 
Ml, 90=0. 408  Model  Aerodynamic  Data  Without  Un- 
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Lift  Coefficient  vs.  Control  Plane  Angle  of  Attack 
Fiberglass  Blade,  0 Twist,  = 5° , h = 0.299, 
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Figure  31. 


Power  Coefficient  vs.  Lift  Coefficient  for 
Collective  Pitch  Sweep,  Fiberglass  Blade, 

0 Twist,  = 5^^  ^ = 0.299,  gQ  = 0.408. 
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Figure  35.  Rotor  Drag  vs.  Collective  Pitch  for  Full-Scale 


Rotor,  -8°  Twist,  u = 0.300, 
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Figure  36. 


Rotor  Power  vs. 
Rotor,  -8°  Twist 


Collective  Pitch  for  Full-Seal 
, H = 0.300,  M,  on  = 0.740. 
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Figure  37. 


Lift  Coefficient,  C /a  x 10 

Power  Coefficient  vs.  Lift  Coefficient  for 
Full-Scale  Rotor,  -8*’  Twist,  u --  0.300, 

Ml, go  = 0.740. 
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Figure  38.  Experimental  Lift  Coefficient  vs.  Calculated 
Coefficient  for  Full-Scale  Rotor,  -8°  Twist, 
p = 0.300,  gQ  = 0.740. 
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-8  Twist,  M = 0.300,  M^^gg  = 0.740.  ® 
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Figure  50. 


Experimental  Power  Coefficient  vs.  Calculated 
Power  Coefficient  Using  Modified  Airfoil 
Tables,  Fiberglass  Blade,  -8°  Twist. 
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Figure  52.  Measured  and  Calculated  Chord  Bending  Moment  Time 
Histories,  Fiberglass  Blade,  0 Twist 

= 0.399,  M,  an  = 0-^34,  a = 0.5°  (Cond.  25). 
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Figure  57.  Measured  and  Calculated  Beam  Bending  Moment  Time 
Histories,  Fiberglass  Blade,  -8°  Twist, 

M = 0.399,  M^  gQ  = 0.434,  = 0.5  (Cond.  25). 
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Figure  63.  Measured  and  Calculated  Beam  Bending  Moment  Time 
Histories,  Fiberglass  Blade,  0°  Twist,  = 5 
la  = 0.399,  M,  qh  = 0.434,  a = 0.5°  (Cond.  25). 
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Figure  75.  Measured  and  Calculated  Beam  Bending  Moment  Harmonics,  Fitergla 
Blade,  0°  Twist,  M = 0.399,  M,  qq  = 0.434,  = 0.5  (Cond.  25) 
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Figure  76.  Measured  and  Calculated  Chord  Bending  Moment  Harmonic 
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Figure  80.  Measured  and  Calculated  Chord  Bending  Moment  Harmonics,  Fiberglass 
Blade,  0°  Twist,  m = 0.299,  M-,  qp  = 0.408,  cx  = 0°  (Cond.  68). 
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Measured  and  Calculated  Chord  Bending  Moment  Harmonics,  Aluminum 
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Figure  101.  Continued 


r/f) 


r/R  = 0.20 
A--  Measured 

Computed  Without 
Unsteady  and  Yawed 
Flow  Effects 
0--Computed  With 

Unsteady  and  Yawed 
Flow  Effects 


0 45  90  135  180  225  270  315  360  405  450 

Azimuth,  ieg 


Figure  103.  Example  oi  Effect  of  Unsteady  Aerodynamics  an: 
Yawed  Flow  Effects  on  Chord  Bending  Moment 
Fiberglass  Blade 


^ - Measured 
♦ --Computed  Without 
Torsion 

0"“ Computed  With  Torsion 


0 45  90  135  180  225  270  315  360  405  450 

Azimuth,  deg 

Figure  104.  Example  of  Effect  of  Addition  of  Torsional  Mod 
Shape  Fiberglass  Blade,  -8°  Twist,  p = 0.502. 
^1,90  “ 0.467,  = 0.5  (Cond.  25). 
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Chord  Bending  Moment,  in. -lb 


r/R  = 0.20 
A--Measured 
■»-- Computed  Without 
Inplane  Spring 
0--Computed  With 
Inplane  Spring 


Azimuth,  deg 


Example  of  Effect  of  Inplane  : 
Blade,  0°  Twist,  p = 0.399,  M. 

a = 0.5°  (Cond.  25). 


Beam  Bending  Moment 


. M^_go  . 0.43/,  a,^=’o.5“  (Cond.  25). 


Azimuth,  deg 


Figure  107.  Three  Successive  Revolutions  of  Chord  Bending^ 

Moment  Time  History,  Fiberglass  Blade,  -8  Twist 
!j  = 0.399,  M,  an  = 0.434,  a = 0.5®  (Cond.  25). 
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TABLE  1.  H-34  FIBERGLASS  MODEL  ROTOR  BLADE 

STIFFNESS  PROPERTIES 


Ar 

(in.  ) 


in 

(in. ) 


I 


(in%) 


(in!^) 


J 

(in.^) 


.23 

3.00 

.466 

2.68 

2. 

37 

.52 

3.23 

.578 

3.54 

.28 

3.75 

1.870 

5.32 

.30 

4.03 

.268 

.268 

.45 

4.33 

1.600 

1.600 

2 

37 

.62 

4.78 

1.030 

1.030 

1 

36 

.63 

5.40 

.814 

.814 

1 

36 

.22 

6.03 

1.120 

1,129 

1 

99 

„28 

6.25 

.330 

.359 

52 

.54 

6.53 

.336 

.332 

50 

.56 

7.07 

.1907 

.1497 

36 

.65 

7.63 

.0711 

.0712 

0440 

.125 

8.28 

.003148 

.02190 

0070 

.495 

8.405 

.00916 

.01730 

00427 

3.630 

8.900 

.001236 

.01240 

00370 

.270 

12.53 

.001236 

.01240 

00343 

41.035 

12.80 

.001060 

.01152 

00319 

.320 

53.835 

.00057 

.00625 

00174 

.845 

54.155 

.00057 

.00625 

00174 

NOTE:  The  following  values  of  E and  G should  be  used  to 

obtain  bending  and  tv/isting  stiffnesses  for  the 
fiberglass  blades: 


E = 2.50  X 10^  Ib/in.^ 
G = 1.01  X 10^  Ib/in.^ 


i 


j 

1 


.1 

I 
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TABLE  2. 

H-34  ALUMINUM  MODEL 
STIFFNESS  PROPERTIES 

ROTOR  BLADE 

hr 

(in. ) 

^in 
(in. ) 

^x 

(in.^) 

I 

(1^4) 

J 

(in 

.23 

3.00 

.117 

.670 

.600 

.52 

3.23 

.145 

.885 

.28 

3.75 

.468 

1.330 

.30 

4.03 

.0670 

.0670 

.45 

4.33 

.400 

.400 

.600 

.62 

4.78 

.258 

.258 

.344 

.63 

5.40 

.204 

.204 

.344 

.22 

6.05 

.280 

.280 

.503 

.28 

6.25 

.0825 

.0898 

.132 

.54 

6.53 

.0840 

.0830 

.127 

.56 

7.07 

.0477 

.0374 

.0911 

.65 

7.63 

.05456 

.0574 

.0367 

.125 

8.28 

.002415 

.0177 

.00583 

.495 

8.405 

.001504 

.,0140 

.00356 

3.630 

8.90 

.000948 

.0100 

.00308 

.270 

12.53 

.000948 

.0100 

.00286 

41.035 

12.80 

.000813 

.00929 

.00266 

.320 

53.835 

.000437 

.005041 

.00145 

.845 

54.155 

.000437 

.005041 

.00145 

NOTE : The 

following  values  of  E and  G should  be  used 

to  obtain 

bending  and  twisting  stiffnesses  for  the  aluminum  blade: 
E = 10.0  X 10^  Ib/in.^ 

G = 4.00  X 106  Ib/in.^ 

The  values  of  I . I . and  J between  r = 8.28  in.  and  r = 
X y 

53.835  in.  have  been  corrected  to  reflect  the  stiffness 
test  results  for  the  aluminum  blade. 


TABLE  3. 


H-34  MODEL  kOl'OK  BLaDE  1.,,EKTIA 
AMD  CENTER  OF  GRAVITY  DATA 


\ 

Ar 

(in. ) 

^in 
(in. ) 

w 

(Ib/in. ) 

P 

(10‘^ 

I 

X 9 

lb-sec'^) 

Pi 

/ -3,  y 
(lO  ■^Ib-sec'^) 

Y 

(in. 

.23 

3.00 

.147 

.0036 

.177 

33 

[ 

.52 

3.23 

.147 

33 

.28 

3.75 

.475 

10 

.20 

4.03 

.149 

33 

. 10 

4.23 

.162 

.0036 

.177 

.30 

\ 

f. 

.17 

4.33 

.180 

.0801 

. 254 

27 

[■ 

.14 

4.50 

.273 

.0801 

. 254 

18 

\ 

t 

.14 

4. 64 

.343 

.0801 

.254 

14 

i: 

r 

.62 

4.  78 

.343 

.0928 

.266 

14 

.17 

5.40 

.343 

.0892 

.0892 

0 

.0 

i- 

.07 

5.57 

.295 

.27 

J . 04 

.112 

. 12 

5.81 

.169 

.0892 

.0892 

i 

.22 

6.03 

.169 

.0645 

.0645 

>%. 

» 

. -iO 

6.25 

.048 

.0177 

.0177 

1 

.54 

6.53 

.0822 

.0193 

.0199 

.56 

7.07 

.0588 

.0111 

.00955 

|:* 

.65 

7.63 

.0414 

.0042 

.00504 

f' 

.125 

8.28 

.0176 

.00037 

.00314 

■ 

.495 

8.405 

.0081 

.00025 

.00268 

3.630 

8.90 

.0094 

.00013 

.0086 

K.- 

i-; 

.270 

12.53 

.0106 

.00014 

.0089 

1: 

40.135 

12.80 

.0121 

.00018 

.0114 

I 

.72 

52.935 

.0202 

.00034 

.0221 

r 

.18 

53.655 

.0178 

.00034 

.0221 

f 

.32 

53.835 

.0276 

.00034 

.0221 

r 

.845 

54.155 

.0044 

.00009 

.00571 

0 

.0 

NOTES:  1. 


2. 


3. 


Mass  moments  of  Inertia  are  wi  th  respect  to  the 
local  chord  line  and  an  axis  normal  to  it  at  the 
feathering  axis. 

Mass  moments  of  inertia  include  parts  at  the  root 
that  flap,  lag,  and  pitch  with  the  blade.  Parts 
that  flap  and  lag  but  do  not  pitch  are  excluded. 
The  flapwise  position  of  the  blade  center  of  grav- 
ity is  on  the  chord  line  for  all  blade  stations. 


Item 

1 

2 

3 

5 

6 


7 

8 

9 

10 


TABLE  4.  H-34  MODEL  ROTOR  BLADE 

MISCELLANEOUS  DATA 


Parame'cer  or  Description 


Value  Units 


Radius 

Blade  Airfoil  Chord  Outboard  of  8.9" 
Radial  Station 

Blade  Shank  Chord  Inboerd  of  8.9"  Radial 
Station 

Blade  Linear  Twist  From  Center  of 
Rotation  to  Tip 

Blade  Sets  a,  c,  d 
Blade  Set  b 

Blade  Shear  Center  Locations 

Beamwise  With  Respect  to  Feathering 
Axis 

Chordwise  With  Respect  to  Feathering 
Axis 

Location  of  Blade  Hinges 

Radial  Station  of  Coincident  Flapping 
and  Lead-Lag  Hinge 

Feathering  Bearing  Outboard  End  Radial 
Station 

Feathering  Bearing  Inboard  End  Radial 
Stations 

(The  feathering,  flapping  and  lag 
bearing  axes  intersect  at  the  same 
poin. ) 

Location  of  Lead-Lag  Damper 

Rotary  Damper  on  Blade  Lead-Lag  Hinge 

Pitch  Flap  Coupling  Ratio 

Damping  Coefficient  of  Lead-Lag  Damper 
(Expressed  as  a damping  moment 
about  the  hinge) 


55.0 

in.. 

2.69 

in. 

1.19 

in. 

0.0 

deg 

-8.0 

deg 

0.0 

in. 

0.0 

in. 

3.0 

in. 

5.60 

in. 

4.64 

in. 

0.0 

17.0 


Tabs  or  Flaps 

Blade  Set  c Deflection  Over  the  Aft  20%  of 
the  Airfoil  Chord  5.0 


in.  ■ 

lb- 

sec 


deg 
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TABLE  4.  H-34  MODEL  ROTOR  BLADE 

MISCELLANEOUS  DATA  (Cont.'d.) 


Item 


Parameter  or  Description Value  Units 


11  Pitch  Control  Geometry 

At  0 Collective  and  Cyclic  Input: 

Distance  Aft  Along  Flapping  Hinge  Axis 

to  Pushrod  Upper  End 

Radius  on  Swashplate  of  Pushrod  Lower 

End 

Distance  of  Plane  of  Swashplate  Below 
Plane  of  Flapping  Hinges 
Angie  on  Swashplate  Between  Pushrod 
Lower  End  and  Plane  of  Shaft  and  Lead- 
Lag  Axis 

Pitch  Com  -ol  Spring  Rate 
Airfoil  Section  - NACA  0012 


1.43 

in. 

1.95 

in. 

6.31 

in. 

40.5 

deg 

21,000 

in. 

lb/ 

rad 

TABLE  5.  H-34  FIBERGLASS  MODEL  ROTOR  BLADE 

EQUAL  SEGMENT  STIFFNESS  PROPERTIES 


9 

10 
11 
12 
13 
!/■ 

15 

16 

17 

18 

19 

20 


Station 

No. 

OO-lb-in.2) 

El 

(lO^lb-^in.2) 

GJ 

(lO^lb-in.2) 

1 

1.7 

2.59 

194.4 

2 

1.7 

2.59 

194.4 

3 

.460 

.438 

15.6 

4 

.00345 

.0337 

.396 

5 

.00292 

.0302 

.351 

6 

.00265 

.0288 

.322 

,00265 

,00194 


.0288 

.0212 


.322 

.238 


TABLE  6.  H-34  ALUMINUM  MODEL  ROTOR  BLADE 

EQUAL  SEGMENT  STIFFNESS  PROPERTIES 


Station 

No. 

El 

(lO^lb-in.2) 

El 

(lO^lb-in.2) 

GJ 

(lO^lb-in 

1 

1.70 

2.59 

196.3 

2 

1.70 

2.59 

196.3 

3 

.757 

.709 

35.0 

4 

.0106 

.109 

1.30 

5 

,00896 

.0975 

1.67 

6 

.00813 

.0930 

1.07 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 


.00813 

,00595 


.0930 

.0685 


1.07 

.793 


.261 

.0904 

.209 

0 

.0 

.261 

.0904 

.209 

-0 

.178 

.0905 

.0311 

.0311 

0 

.0 

.00989 

.000207 

.00725 

.0105 

.000148 

.00960 

.0121 

.000180 

.0114 

,0121 

,0140 


.000180 

.000223 


.0114 

.0144 


TABLE  7.  H-34  MODEL  ROTOR  BLADE 

EQUAL  SEGMENT  INERTIA  AND 
CENTER  OF  GRAVITY  DATA 


Station 

No. 


-3, 


(Ib/in)  (lO  ■^Ib-sec^) 


(I0"^lb4ec2)  (iSf) 


TABLE  8.  H-34  MODEL  ROTOR  BLADE  MODE  TYPES, 

FREQUENCIES,  AND  INERTIAS 

Fiberglass  Blade,  0°  Twist 


Mode 

No. 


Mode 

Type 


Natural  Frequency 
(per  rev) 


Genera] ized  Inertia 
(lb-ft2) 


1 

B 

1.045 

0.00660 

2 

C 

0.305 

0.00664 

B 

2.588 

0.00455 

4 

C 

3.171 

0.00554 

S 

B 

4.381 

0.00546 

6 

C 

7.538 

0.00614 

Fiberglass  Blade 

, -8°  Twist 

1 

B 

1.045 

0.00660 

2 

C 

0.305 

0.00664 

3 

B 

2.594 

0.00437 

4 

C 

3.164 

0.00535 

5 

B 

4.387 

0.00536 

6 

C 

7.556 

0.00589 

Aluminum  Blade, 

0°  Twist 

1 

B 

1.045 

0.00661 

9 

C 

0.305 

0.00668 

3 

B 

2.769 

0.00495 

4 

C 

4.497 

0.00561 

5 

B 

5.282 

0.00593 

6 

B 

9.000 

0.00606 

B = 

Beamwise  Mode 

Rotor  Speed 

=730  RPM 

C = 

Chordwise 

Mode  Root 

Collective  Pitch  = 12° 

TABLE  9.  CORRELATION  CRITERIA 


A. 

Mean  Load  and  Rotor  Performance 

Excellent 

0/5 

Good 

5/10 

Fair 

10/15 

Poor 

15/ 

B. 

Overall  Amplitude  or 

Phase 

Excellent 

0/10 

Good 

10/20 

Fair 

20/30 

Poor 

30/ 

C. 

Harmonic  Amplitude 

Scale /Harmonic  1 

2 3 

4 

5 

Excellent  0/10 

0/20  0/30 

0/35 

0/40 

Good  10/20 

20/40  30/60 

35/70 

40/80 

Fair  20/30 

40/60  60/90 

70/105 

80/120 

Poor  30/ 

60/  90/ 

105/ 

120/ 

D. 

Harmonic  Phase 

Scale /Harmonic 

1,2 

3,4,5 

Excellent 

0/10 

0/20 

Good 

10/20 

20/40 

Fair 

20/30 

40/60 

Poor 

30/ 

60/ 

TABLE  10.  PRESENTATION  AND  GRADING  OF 

POST  TEST  PERFORMANCE  PARAMETERS 


TABLE  10.  Continued 
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TABLE  14.  MEASUREP  BEAM  AND  CHORD  BENDIN5  MOMENTS, 
FIBERGLASS  BLADE,  0°  TWIST 

Lt  = 0.502,  Ml  an  = 0.467,  = 5.0®  (GOND.  44) 
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TABLE  15.  CALCULATED  BEAM  AND  CHORD  BENDING  MOMENTS, 
FIBERGLASS  BLADE,  0°  TWIST 

M = 0.502,  Ml  9Q  = 0.467,  = 5.0°  (COND.  44) 


TABLE  16.  MEASURED  BEAM  AND  CHORD  BENDING  MOMENTS, 
FIBERGLASS  BLADE,  0°  TWIST 
|a  = 0.299,  Ml  90  = 0.408,  = 0°  (COND.  68) 


SIN  0.0  20.92 

CCS  -0.26  -5.74 


TABLE  17.  CALCULATED  BEAM  AND  CHORD  BENDING  MOMENTS. 
FIBERGLASS  BLADE,  0®  TWIST 

}a  = 0.299,  Ml  90  = 0.408,  = 0®  (COND.  68). 
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TABLE  19.  CALCULATED  BEAM  AND  CHORD  BENDING  MOMENTS, 
FIBERGLASS  BLADE,  -8°  TWIST 

M = 0.399,  M,  = 0.434,  = 0.5  (COND.  25) 
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TABLE  23.  CALCULATED  BEAM  AND  CHORD  BENDING  MOMENTS, 
FIBERGLASS  BLADE,  -8°  TJIST 
|a  = 0.299,  M,  qa  = 0.408,  a = 0°  (COND.  68) 
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TABLE  26.  ^EASURED  BEAM  AND  CHORD  BENDING  MOMENTS 
FIBERGLASS  BLADE,  0°  TWIST,  <5^  = 5° 

M = 0.502,  M]^^90  = 0.467,  = 5.0°  (CO 
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TABLE  28.  MEASURED  BEAM  AND  CHORD  BENDING  MOMENTS, 
FIBERGLASS  BLADE,  0°  TWIST,  6^  = 5° 

M = 0.299,  Ml  90  = 0.408,  a = 0°  (GOND.  68) 
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TABLE  30.  MEASURED  BEAM  AND  CHORD  BENDING  MOMENTS, 
ALUMINUM  BLADE,  0°  TWIST 

M = 0.399,  Ml  90  = 0.434,  £%  = 0.5°  (COND.  25) 


— < O 'J  rT)  r-<  ^ o fM  (<1 


0 o 

1 I 


J U O 'J 


J O J — • 
I I I 


'T>  00  o rn  sj-  00  i»"\ 
■w  n r~  <-r  '-H  m 


j -o 


O r J —<  0 O ■-< 


■J  o 

I I 


(NjO'int^O(\jLnir>r^f«- 
vO  r^O'NinrNj^tor-oo' 

• ••••••••• 

'-^00-«000000 
I I I I I I I 


0 O O O 

1 I I 


0 o o 

1 I 


o o 
I 


lA 


Osl 

00 

>0 

f\l 

>»■ 

sO 

•o 

O' 

o 

o 00 

00 

O' 

O' 

cr 

in 

o 

o 

00 

m 

00 

00 

m 

o 

og 

00 

>0 

rg 

m 

(M  n 

00 

OO 

eg 

rg 

rg 

rg 

in 

rO 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• ■ 

• 

• 

• 

• 

• 

• 

• 

• 

O 

o 

o 

Q 

o 

o 

o 

o 

o 

(\l  r\J 

•o 

-o 

in 

eg 

I I I 


— « 

u. 

LL 

30 

oa 

_) 

-1 

o 

O' 

m 

o 

cr 

U1 

CJ 

-r 

1 

a^ 

rn 

r- 

h- 

U1 

g" 

g- 

O' 

rg 

1 

o 

o 

e- 

(A 

o 

o 

in 

z 

'T 

O 

CO 

rg 

O' 

CO 

4* 

r*H 

r- 

O' 

• 

* 

• 

« 

• 

. 

• 

• 

• 

• 

• 

• 

• 

m 

• 

• 

• 

• 

• 

• 

o 

o 

o 

r~i 

•—4 

•— < 

o 

o 

w 

in 

rg 

o 

rg 

in 

rn 

r'H 

r-g 

1 

1 

1 

1 

1 

1 

T 

1 

>— 

1 

1 

1 

1 

1 

h- 

z 

Z 

UJ 

LU 

X 

z 

o 

o 

eg 

e- 

O' 

in 

m 

r~ 

CO 

jO 

f^ 

o 

X 

GO 

LA 

O' 

o 

CO 

rg 

e- 

OO 

■o 

X 

m 

go 

"4- 

00 

00 

O' 

«— 1 

g 

00 

TO 

en 

o in 

00  'S’ 

e» 

o 

in 

sO 

eg 

• 

• 

• 

• 

« 

• 

• 

• 

< 

• 

Q 

t 

• 

• 

• 

t 

• 

• 

* 

X 

in 

en 

OO 

in 

00 

00 

eg 

go 

oc 

n 

. 

eg 

o 

•—4 

o 

m 

<t 

1 

1 

1 

1 

1 

o 

1 

1 

1 

r-H 

1 

1 

fH 

OJ 

X 

1 

1 

1 

3) 

u 

Q 

Q 

UJ 

O' 

•H 

eg 

00 

rg 

-0 

en 

LU 

m 

m 

■a 

m 

r* 

o 

o 

h- 

o 

QC 

CM 

o 

m 

eg 

GO 

n- 

f'» 

30 

QC 

eg 

eg  O' 

'4- 

m 

't 

4“ 

in 

O 

eg 

X 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

X 

• 

• 

t 

• 

e 

• 

• 

• 

• 

• 

o 

'0 

o 

o 

rn 

o 

eg 

o 

in 

0 

2 

O 

m 

o 

tn 

't 

m 

m 

< 

1 

1 

1 

1 

1 

1 

< 

1 

1 

OJ 

lU 

X 

X 

rg 

O' 

■4- 

rg 

m 

r- 

eo 

o 

'<■ 

eg 

O >}• 

O' 

rg 

in 

rg 

in 

h- 

fH 

GO 

g- 

o 

eg 

eg 

g 

fH 

-t  e- 

eg 

m 

r^ 

f*4 

4- 

vO 

•O 

!^ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• * 

• 

t 

• 

• 

• 

• 

• 

• 

eg 

F^ 

4- 

o 

in 

O 

4* 

F^ 

eg 

eg  CO 

vO 

eg 

m 

>r 

GO 

'T 

go 

in 

I -<  I -*  I rg  I 


TC 

a: 

< 

X 


to 


gj 

>r 

n 

g) 

rg 

n 

g- 

00 

o 

n 

o 

eg 

a 

eg 

o 

f-H 

o 

go 

o 

o 

in 

o 

g- 

o 

eg 

o 

fH 

o 

m 

• 

• 

• 

• 

• 

• 

• 

• 

« 

• 

t 

• 

• 

• 

• 

• 

• 

■ 

• 

• 

o 

eg 

1 

n 

eg 

o 

e- 

o 

00 

o 

30 

o 

f\j 

r> 

g- 

o 

in 

n 

go 

o 

in 

2 

m 

3 

z 

m 

z 

1/) 

z 

00 

X 

z to 

z 

00 

z 

00 

z 

00 

z 

00 

o 

M4 

o 

o 

»— < 

a 

>—* 

o 

QC 

l-H 

o 

O 

O 

o 

FH 

a 

i/) 

o 

1/) 

u 

in 

o 

1/} 

o 

o 

< 

1/) 

u 

1/) 

O 

l/> 

i_> 

to 

o 

to 

o 

W1  N1 


< 

(- 

in 


W1  tf«  M 

in  m m 

'Oj  'Tj 


249 


•aw;ff2i25i£6 


• 

r-H 

O' 

nT 

o 

o 

0 

f" 

fNJ 

f H 

in 

in 

3D 

♦t 

o 

03 

-o 

o 

-o 

r-H 

Ln 

r^ 

rH 

o 

in 

•o 

o 

o 

O 

o 

o 

r— 4 

in 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

(N 

o 

o 

o 

o 

O 

o 

O 

o 

•D 

o 

O 

o 

"D 

o 

D 

O 

a 

o 

u 

o 

DO 

o 

00 

in 

'i- 

rg 

4- 

4 

m 

in 

i-H 

•o 

3D 

-o 

00 

4- 

O' 

00 

u^ 

vj- 

r- 

rg 

fO 

-0 

O' 

rr> 

C* 

u~\ 

■o 

4 

nj 

•o 

o 

lA 

-0 

r\J 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

0 

• 

• 

ft 

f 

ft 

0 

f\j 

(M 

o 

o 

»— 4 

rO 

O 

o 

rH 

'■3 

r— H 

O 

J 

3 

2 

-J 

rH 

\l 

o 

O' 

l-J 

«n 

in 

o 

-n 

u^ 

1 

*NJ 

tn 

0^ 

g- 

-~H 

3 

D 

ig 

w e 

1 

4- 

■O 

J 

rg 

'J 

o 

4 

r-^ 

O' 

U 

4 

nT 

4 

UD 

O' 

0- 

O' 

4 

u^ 

.3 

CQ  d 

ft 

4 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

a 

o 

o 

f~H 

<— j 

w 

0 

r*J 

Q 

r\l 

1-^ 

n 

.3 

rn 

0 

CJ 

Q 

■— 

1 

1 

1 

1 

K- 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Pi  H <t 

1— 

2. 

O c/D  CO 

2 

lU 

X M <f 

l-J 

757 

CJ  s . 

i. 

3 

^ o 

O 

■o 

t 

O' 

•0 

4 

00 

'll 

OJ 

JD 

•o 

X 

rg 

fNJ 

N. 

m 

7N 

CO 

g- 

rg 

0 

<n 

Q 

2: 

ro 

4 

(M 

O' 

0 

4 

O' 

4 

m 

O' 

CO 

m 

0 

0 

r— J 

m 

3 

0 

4 

2.C-  II 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

0 

ft 

ft 

ft 

ft 

ft 

• 

ft 

ft 

ft 

ft 

<■  o 

2i: 

m 

o 

m 

CD 

CJ 

fO 

OJ 

ro 

r-J 

PC 

*—4 

IM 

CJ 

4 

CO 

m 

in 

IT 

CJ 

3 

o 

< 

1 

0 

! 

1 

1 

1 

1 

•»  (TN 

OJ 

X 

< u 

cO 

3 

W Q 

PQ  < 2 

a 

Q 

_) 

■u 

4 

in 

-0 

CJ 

lA 

ro 

O' 

O 

■o 

ai 

O' 

0 

OD 

ro 

N- 

0* 

O' 

r-4 

in 

m 

Q PQ  ^ 

h- 

rvj 

D 

O 

fO 

in 

O' 

J5 

f-4 

4 

eg 

4 

J— 

CJ 

4 

0 

0 

eg 

4 

in 

m 

0 

O' 

m 

W CA 

3 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

ft 

3 

ft 

ft 

ft 

ft 

• 

ft 

ft 

ft 

• 

ft 

H 2 

o. 

CJ 

rn 

(O 

CJ 

CJ 

r-H 

CJ 

CJ 

»— H 

Cl 

Q. 

0 

0 

^■4 

0 

0 

0 

3 

0 

< S CO 

2" 

1 

1 

X 

1 

1 

1 

1 

1 

1 

r-I2  . 

3 

3 

:d  M o 

•3 

13 

O 2 

r-H 

o 

3 

ro 

4 

in 

r'^ 

0 

-n 

c- 

0 

O' 

r“H 

in 

u\ 

m 

<j  < a 

r— 4 

o 

■o 

in 

O' 

CO 

O 

CO 

o 

rO 

OJ 

4 

3 

'T' 

»n 

0 

i—J 

O' 

4 

■ :> 

3 

o 

m 

O 

rn 

Cl 

■ 4 

m 

rg 

0 

m 

4 

TABLE  32.  MEASURED  BEAM  AND  CHORD  BENDING  MOMENTS, 
ALUMINUM  BLADE , 0°  TWIST 

|i  = 0.502,  Mi^9o  = 0.467,  am  = 5.0°  (COND.  44) 
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TABLE  34.  MEASURED  BEAM  AND  CHORD  BENDING  MOMENTS, 
ALUMINUM  BLADE,  0°  TWIST 

H = 0.299,  M,  pp,  = 0.408,  a = 0°  (COND.  68) 
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LIST  OF  SYI-IBOLS 


n 


a 


Bn 

b 

''D 

Co/a 


Cp/o 

Cp^/o 

CPp/o 

Cq/o 

Cj^/O 

Cy/O 

E 

I 

^EQ 

1 


L 

M 

M 


Harmonic  coefficient  of  cosine  terms 

Vertical  intercept  of  regression  line 

Harmonic  coefficient  of  sine  terms 

Slope  of  regression  line 

Drag  coefficient 

Rotor  drag  coefficient 

Lift  coefficient 

Rolling  moment  coefficient 

Maximum  lift  coefficient 

Pitching  moment  coefficient 

Yaw  moment  coefficient 

Rotor  power  coefficient 

Induced  power  coefficient 

Parasite  power  coefficient 

Rotor  torque  coefficient 

Rotor  lift  coefficient 

Rotor  side  force  coefficient 

2 

Modulus  of  elasticity,  Ib/in. 

4 

Area  moment  of  inertia,  in. 

Equivalent  area  moment  of  inertia,  in.^ 
Segment  length,  in. 

Rotor  lift,  lb 

Bending  moment  magnitude,  in. -lb 
Mach  number 


255 


Advancing  tip  Mach  number 
Quality  factor 
Blade  radius,  in. 

Multiple  correlation  coefficient 
Blade  radial  location,  in. 

Reynolds  number 
time,  sec 

Blade  segment  weight  distribution,  Ib/in. 

Distance  from  blade  segment  quarter  chord  to  blade 
segment  center  of  gravity,  in.  (positive  forward) 

Angle  of  attack,  deg 

Mast  angle  of  attack,  deg 

Tip  path  plane  angle  of  attack,  deg 

Flap  angle  on  trailing  edge  flap  airfoil,  deg 

Advance  ratio 

Blade  segment  mass  moment  of  inertia  in  out-of- 
plane direction,  lb-sec2 

Blade  segment  mass  moment  of  inertia  in  in-plane 
direction,  lb-sec^ 
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